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Abstract 

The complex-shift method is apphed to the Kuzmin-Toomre family 
of discs to generate a family of non-axisymmetric flat distributions of 
matter. These are then superposed to construct non-axisymmetric flat 
rings. We also consider triaxial potential-density pairs obtained from 
these non-axisymmetric flat systems by means of suitable transforma- 
tions. The use of the imaginary part of complex-shifted potential- 
density pairs is also discussed. 
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1 Introduction 

One of the difficulties encountered in modelling self-gravitating systems is set 
by the potential theory. In general, to calculate the gravitational potential 
associated with a matter distribution of arbitrary shape, one has to use 
numerical methods. On the other hand, it is convenient to have a set of 
analytical potential-density pairs that provide a description of gravitating 
systems, in particular, for those that deviate from spherical symmetry. A 
survey of such models can be found in [l], chapter 2. 

A simple way to construct axisymmetric potential-density pairs that rep- 
resent discs is the image method that is usually used to solve problems in 
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electrostatics. This method was employed in the context of the Newtonian 
gravity first by Kuzmin [2] , who derived a very simple potential-denstity pair 
of a disc (Section [3]). Several years later, Evans & de Zeeuw [3] showed that 
any axisymmetric disc can be constructed by superposing simple Kuzmin 
discs with different weights. The image method has also been adapted to and 
used in General Relativity to generate several exact solutions of Einstein's 
equations that represent discs (see e.g. [H El El El El EJ [lOl flT]). 

Another technique to construct new solutions in the Newtonian gravity 
was introduced by Appell [l2l [13] who considered the complexification of the 
potential of a point mass. This method, known as complex-shift method, was 
also used in General Relativity [Ml [HI [16] , but see some misinterpretations 
in the Appendix of [17j. More recently, the complex-shift method has been 
applied to spherical and axisymmetric Newtonian parent systems to generate 
new analytical axisymmetric and triaxial potential-density pairs [181 US] • 

In this work, we use the complex-shift method to derive potential-dentity 
pairs for non-axisymmetric flat discs and rings and their "inflated" versions. 
In Section [21 we summarize the method of complexification. In Section [3l we 
complex shift the Kuzmin- Toomre family of discs to derive a family of non- 
axisymmetric flat discs. The members of this class are then superposed to 
generate further non-axisymmetric flat distributions of matter that can be 
interpreted as complex-shifted rings. This is done in Section [H In Section [5l 
we use a general transformation to "inflate" the previous non-axisymmetric 
flat systems to generate triaxial potential-density pairs. We derive a gen- 
eralized pair that includes an example studied by Ciotti & Marinacci [19] . 
In Section [6l we show that the imaginary part of a complex-shifted system 
can be superposed to another potential-density pair to generate physically 
acceptable systems which are not symmetrical with respect to a coordinate 
axis. At last, in Section [7] we summarize our results. 

2 The complex-shift method 

In this section, we briefly expose the idea of the complex-shift method. The 
formalism is similar to the presented by Ciotti & Giampieri [18j and Ciotti & 
Marinacci [19]. We start by considering a gravitational potential '&(x) and 
an associated density distribution p(x), where x = {x,y,z) is the position 
vector. We also assume that p(x) satisfies the Poisson equation 

= AnGp. (1) 
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The basic idea of the complex-shift method is to obtain new potential-density 
pairs by a displacement of a given potential-density pair (the parent system) 
on the imaginary axis. The complexified potential is defined as $c = '^'(x — 
ia), where i is the imaginary unit and a is a real shift vector. Since the 
complex shift is a linear transformation and the Poisson equation is a linear 
partial differential equation (PDE) acting on the x vector, it follows that 

V'$e = ^TrGpc, (2) 

where pc = p(x — ia) is the shifted mass density. By separating the real and 
imaginary parts of and pc, one obtains two real potential-density pairs. 

Unfortunately, there is no guarantee that the resulting pairs describe 
physically acceptable (i.e. everywhere positive) gravitating systems. In fact, 
Ciotti & Giampieri [18] have shown that the imaginary part of the complex- 
shifted density always changes sign because the total mass of the complexi- 
fied potential-density pair coincides with the total mass of the seed density 
distribution, and thus / 9(/3c)d^x = 0. However, the real part can be posi- 
tive everywhere. This usually happens for a finite domain of the shift vector 
a. 



3 Thin shifted Kuzmin-Toomre discs 

In this Section, we apply the complex shift to the Kuzmin-Toomre family 
of discs O [20] to generate flat distributions of matter without axial sym- 
metry. We begin with the first member of this family, which was the simple 
model derived by Kuzmin. The potential in the usual cylindrical coordinates 
(i?, z, ip) is given by 

GM 

^ = =, 3) 

where 6 is a non-negative parameter. Using Poisson equation ([1]), the dis- 
continuous normal derivative on z = introduces a surface density of mass 



"^^^ = 2^ 97 



(4) 

2=0 



which for the potential ([3]) results 



. . bM 

a{R) = ry^. 5 

27r {m + 62)3/2 

The result of a complex shift of an axisymmetric potential depends not 
only on the length but also on the direction of the shift vector a. We are 
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interested in the equatorial shift of the potential ([3]), where without loss of 
generality we make a shift along the x-axis, a = (a, 0,0). The complexified 
potential ([3|) reads 

$c = , . (6) 

V(x-m)2 + y2 + (5+|^|)2 

To evaluate the square root, we denote cP = + — + (6+ — 2iax = 
ite*^. Thus, 

U = \df = \j [x2 + y2 _ ^2 + (5 + |^|)2]2 ^ 4^23.2^ (7) 

x'^ -a^ + {h+\z\f . 2ax 

cos 8 = , smt^ = . (8) 

u u 

For simplicity, we restrict to the case cos > everywhere, so we must have 
a < b. This restriction will be assumed henceforth. To make the square root 
a single valued function, we choose d = y^e*^/^, with 



9 /1 + COS6' e sine 

With the help of equations ©-([O]), the real and imaginary parts of the 
potential ([6]) can be expressed as 



\/ + ^cfix'^J X + + ^o?x'^ 

where x = a:^^ + — + + The respective surface densities can be 
calculated using relation dH) 



^{cJc) = 1 ^ ^ . 2^ - Ve2 + 4a2x2 , (12) 

2V27r (^2 + 4^23.2)3/2 V J 



Mabx 2i + 

9(^c) = ^ ' = , (13) 

+ 4a2x2)3/2 M + ^^2 + 4„2^2 



where we defined ^ = + + 6^ — a^. Equations (fT2]) and (fT3]) may also 
be obtained from the real and imaginary parts of the complexified surface 
density ([5]). 
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As was stated in Section [2l the imaginary part of the shifted Kuzmin- 
Toomre disc changes sign and becomes negative for x < 0, so we focus on 
the real part. The expansion of the density (jl2p about the origin yields 



Mb 3Mfe(6^ + 4a^)x^ 3Mby^ 

^'^'^ ^ 27r(62 - a2)3/2 " 4vr(62 - a2)7/2 " 4vr(62 - a2)5/2 " ^ ^ 

The isodensity curves are ellipses with major axis on the y-axis, minor-to- 



major squared axis ratio (b'^—a^)/ (b'^+Aa^) and eccentricity e = ^/Sa"^/ (62 -)- 4a2). 
Some curves of constant density in units of M/b"^ are displayed in Figs.dja) 
and (b) as functions of x/b and y/b for (a) a/b = 0.5 and (b) a/b = 0.8. We 
have a monotone decreasing surface density in the interval < a/b < 0.59, 
whereas the density is always non-negative for < a/b < 0.86. Expanding 
the potential ([TO]) about the origin in the plane z = results 



^ GMjb' + 2a2)x2 GMy^ 

-^"^^'^ .JW^ 2{b^-a^fl^ ^ 2(62_a2)3/2' ^'^> 

which is the potential of an anisotropic harmonic oscillator with the ratio of 
frequencies 

w.x /62 + 2a2 



Therefore, we expect that test particles near the centre move in box orbits 
[1]. The contours in Fig. [l][|b) suggest that the simple potential-density pair 
pop and (jl2p may be used to represent thin galactic discs with a central 
bar. 

The potential-density pair obtained form the real part of the shifted 
Kuzmin-Toomre disc is in fact the first member of the family of non-axisymmetric 
thin disc- like structures. The other members can be calculated by using the 
recurrence relations for the Kuzmin-Toomre family of discs [211 [5] 

^c(n+l) = ^dn) - ^^/^^^cin), (17) 

b d , , 

o-c(n+i) = ~ 2n + l 96^^^"^' '^ = 0>l>---> (18) 

where *l*c{o) ^"^^ '^c{o) given by equations (fTUj) and respectively. For 
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Figure 1: (a) and (b) Isodensity curves of the surface density (fT2]) as func- 
tions of x/b and y/h for (a) a/b = 0.5 and (b) a/b = 0.8. (c) and (d) 
Isodensity curves of the surface density (f20]) as function oi x/b and y/b for 
(b) a/h = 0.4 and (b) a/h = 0.55. 
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instance, the potential and density of the next member read 



+b{b+\z\) 2x- Vx2 + 4a2x2 , (19) 



2V27r(e2 + 4a2x2)5/2 



^c{i) = , , . , ( 3^' - 2e + 4a2x2 - 4a^x^) , (20) 



where x = 2:^ + — a2 + (6 + |z|)2 and ,^ = x2 + y2 _|_ ^2 _ ^2^ ggj^^-^g contours 
of the surface density (j20]) are shown in Figs. [T^c) and (d) as functions of 
x/b and y/6 for shift parameters (a) a/b = 0.4 and (b) a/b = 0.55. Now a 
monotone decreasing density is found in the range < a/b < 0.44, and a 
non-negative density in the interval < a/b < 0.58. 

4 Non-axisymmetric thin rings 

We now use the family of non-axisymmetric thin discs discussed in the 
previous section to construct structures that can be considered as non- 
axisymmetric flat rings. The idea is to use superpositions of different mem- 
bers of the Kuzmin-Toomre family of discs to generate rings (zero mass 
density at the origin), as was done recently by Vogt & Letelier [22]. Let us 
briefly recall the main point of that work. The general expression for the 
surface density of the nth-order Kuzmin-Toomre disc is [21j 



(2n + l)62"+iM , , 

^(") = 2vr(i?2 + 62)n+3/2' ^ = 0,1,..., (21) 



For convenience, we take discs with mass 



2TTb^ai , , 

M = 22 

2n + 1 ' ^ ' 



where ai is a constant with dimensions of surface density. Equation (121)1 is 
then rewritten as 



- {R2 + 52)n+3/2 " (1 + ^2/52)^+3/2 ' ^^^^ 
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Now we consider the following superposition: 
m 

= Z^^k (-1) Cr(n+m-k) 



k=0 

?2\ 



A:=0 

where = m!/[(m — Using 

2m 



62 



(24) 



E^n-ir-Mi + ^ = f , (25) 

fc=0 ^ / \ / 

equation (p^ takes the form 

^ (fi/6)'"^ 

(1 + i?2/?>2)"+™+^/2 

We have that ()26p with m = 1, 2, . . . and n = 0, 1, . . . defines a family of flat 
rings (characterized by zero density on i? = 0). 

It would seem that a complexification of the density (j26p would give 
rise to other ring structures. However, consider a shift it(™''"')(x — ia) and 
evaluate the resulting complex density at the origin. We get 

(l-a2/62)«+™+3/2' ^ ^ 

which is non-zero and even negative when m is odd. This indicates that in 
order to generate non-axisymmetric ring structures we should consider the 
superposition of complex-shifted Kuzmin-Toomre discs. We go back and 
apply a shift to equation (p3]) . evaluate it at (x = 0,y = 0), 

= (l_,2j52)-3/2' (28) 

and write the superposition 

£ (1 - a2/62)™-'=+"+3/2 = _ ^2/^2)-+n+3/2 E " ^2 J " ^^^^ 

Now we impose the condition that the sum (j29p be zero. Noting that 



m r / 9\ -\in—k / '?\ k / 9 9\"i 



A;=0 



a' 



'-IP 
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and comparing with (|29|) . we get 



-1 



tm—k 



o \ m—k 
62 



(31) 



If we denote the surface density of the nth-order shifted Kuzmin-Toomre 
disc by crc(„), then the superposition 

2 \ nri—k 



(J. 



(m,n) 



k=0 



^c{n+m—k) 1 



(32) 



win be the equivalent to complex-shifted "rings" with zero surface density 
at the origin, and if a = the superposition reduces to the sum (p^ . The 
corresponding potential is calculated by a sum similar to ([32]) . Note that 
the complex-shift of the family of rings given by (j26p and the superposition 
of complex-shifted Kuzmin-Toomre discs given by (j32p result in different 
structures, because the coefficients used in the sums are different. 
As an example, we analyse the member with m = 1, n = 

,2^ 



(1,0) 



^c(O) 



1 



a 



f^c{l)- 



(33) 



Using equations (fTO|) . (fT2|) and (fT9|) - (pO]) . the explicit expressions for the 
surface density and potential can be cast as 



cr. 



(1,0) 



^/2(e2 + 4a2x2)5/2 



{f + 4a^x^) 2e - Ve + 4a2x2 



3^2 - 2^ + 4a2x2 - 40^x2 



(34) 



$(1,0) 



3V2(x2 + 4a2x2)3/2 



l-^j6(6 + |z|)(2x- Vx' + 4a2. 



(35) 



where and x were defined in Section [3j Near the centre, the approximate 
expressions for the potential-density pair are 



.(1,0) 



cj/63(62 + 9a2)x2 



+ 



$(1,0) 



(?,2_o2)7/2 (52_a2)5/2 



(36) 



7r(T/Gy2 



3(62 _ ^2)1/2 



3(62 _ ^2)5/2 



3(62 -a2)i/2^ 



(37) 
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(a) (b) 



2 




x/b 



Figure 2: (a) and (b) Isodensity curves and surface plot of the surface density 
cTc^'^^/cT/ equation with parameter a/h = 0.5. (c) and (d) Isopotential 
curves and surface plot of the potential ^'i^'^^ / {uiGh) ([55]) with parameter 
a/h = 0.5. 
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which show that the density (potential) has a minimum (maximum) at the 
centre. Figs. [2]^a)-(d) show level curves and surface plots of the dimension- 
less density a^'^^ /ai and potential ^^c''^'' /{aiGb) with shift parameter a/b = 
0.5. Note that we do not have exactly a homogeneous mass concentration 
along a ring. Two points of maximum occur at (x = ±Xs, y = 0), where Xg is 
the non-zero root of the equation da^'^^ /dx = and has to be found numer- 
ically. Two saddle points also occur at (x = 0,y = ^^J2{b'^ — a?)/?)). For 
a/b = 0.5, we find Xs/b « 0.43 and saddle points at y/b = ±^2/2 « ±0.71. 
The potential shows two minimum points at (x = ±Xp,y = 0), where Xp is 
the non-zero root of the equation d^^'^^ / dx = 0, and two saddle points at 
{x = 0,y = ±(62 - a2)/Va2 + 262). For a/b = 0.5, we find Xp/b ^ 0.42 and 
saddle points at y/b = ±1/2. 

The family of potential-density pairs derived in this Section may be use- 
ful to model ring galaxies, or R galaxies [23], objects in the form of approx- 
imate elliptical rings; in particular a subclass of such galaxies, designated 
RE galaxies, which have crisp, elliptical rings with photographically empty 
interiors. 

5 Triaxial potential- density pairs from complex- 
shifted discs 

The systems studied in previous sections were all flat, i.e., had infinitesimal 
thickness. One way to construct three-dimensional potential-density pairs 
from flat gravitating systems is to apply a transformation that involves the 
2;-coordinate. For instance, Miyamoto & Nagai [21] proposed a transforma- 
tion to "inflate" the thin Kuzmin-Toomre discs to obtain three-dimensional 
potential-density pairs for the disc part of galaxies. Other transformation 
based on a class of even polynomials was used by Gonzalez &; Letelier [lOj 
and Vogt & Letelier [11] to generate thick General Relativistic discs. 

We shall denote by h{z) a general even function of the z-coordinate, 
and start with the real part of the potential of the complex-shifted Kuzmin- 
Toomre disc, where we replace |z| by h{z) 



where now we define C = x'^ + y'^ — a"^ + [b + h{z)]'^. The resulting mass 
density follows directly from Poisson equation ([T]) in the usual Cartesian 




(38) 
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coordinates. We obtain 



4V2^(C2 + 4a2x2)3/2 



M(l-/.2,) Vc+x/C^ + 4a2x2 

H p ^75 3(6 + h) 

4^271 (C2 + 4a2x2)5/2 L 

(SC^ - 2C^C + 4a2x2 - 40^x2) - (C^ + 40^x2) (^2C - + 4a2 



(39) 



Note that to satisfy the restriction cos > in (l8|) , we must have C > 
everywhere. In the case of a Miyamoto & Nagai transformation [23] , h{z) = 
\/ ■\- (? with c > 0, and the potential-density pair (f38l) and (f39l) is the 
same as studied by Ciotti & Marinacci [19], who performed an equatorial 
shift of the Miyamoto-Nagai disc. We can also use the class of polynomials 
used in [ini E] 



h{z) 



—z + C, z < —s, 

Az^ + Sz2"+2, -s<z<s, (40) 
^z + C, z > s, 



where 



2n + l-sw _ sw - 1 _ s{2n + l + sw) 

~ ' ~ 4n(n + l)s2"+i ' 4(n + 1) ' ^ ^ 

with n = 1,2,...; s is the half-thickness of the disc and w is the jump of h^zz 
on z = lbs. For \z\ > s, we have h^z = il and ^,^2 = 0, so there is no matter 
in this region. We will consider the special case w = 0, so h^^zi^ = ^s) = 
and the density also vanishes there. Figs. [3]^a)-(c) show isodensity curves 
in the three orthogonal coordinate planes of the mass density p^-^Q^/ {M /h^) 
for the polynomial with n = 1, half-thickness s/h = \ and shift parameter 
a/h = 0.6. We find that the mass density is non-negative in the range 
< a/6 < 0.62, and when the half-thickness is reduced to s/h = 0.5 the 
interval becomes < a/6 < 0.70. The main qualitative difference between 
our system and the one studied by Ciotti & Marinacci [19] for the Miyamoto- 
Nagai shifted disc lies in the finite thickness of the former. 

The "inflated" version of the ring-like non-axisymmetric flat systems (134p 
and ([35]) also can be found by replacing \z\ by h{z) in the potential ([35]) and 
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Figure 4: Isodensity curves in the three orthogonal coordinate planes of the 
surface density p^'^^ /{ai/b) equation (j42p for shift parameter a/b = 0.5 and 
a Miyamoto-Nagai transformation with parameter c/b = 0.5. 
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calculating the resulting mass density with Poisson equation. One ends with 



= — { 36 + h{2 + aV6' C' + 4a\2 

6V2(C2 + 4a2x2)5/2 / vs ; 

2C - VC^ + 4a2x2) - 36(6 + hf{l - j}?) (?>C - 2C VCM^4^ - 40^x2) | 



cT/62 (1 - h\\ a/c + VC^ + 4a2x2 ^ ^ 

X (C^ + Ao?x^) (?>C - + 4a2x2 - 40^x2) - 156(6 + hf{\ - j}?) 

X (^4C^ - 3CVC^ + 4a2x2 + 4a2x VC^ + 4a2x2 - leCa^x^^ 

-(2 + aV^') (C' + 4a2x2)' (2C - TcM^i^) } . (42) 

Some isodensity curves in the three orthogonal coordinate planes of the mass 
density p'c''^^ /{cri/h) are shown in Figs. [ll|a)-(c), for the case of a Miyamoto- 
Nagai transformation h{z) = \/ + c^. The values of the parameters used 
were a/6 = 0.5 and c/6 = 0.5. As expected, one has a non-axisymmetric 
toroidal mass distribution. 



6 Using the imaginary part of complex-shifted discs 

In Section [21 we commented that the imaginary part of a complex-shifted 
density always changes sign [18] , and therefore is not adequate for modelling 
gravitating systems. In fact, the imaginary part of the shifted Kuzmin- 
Toomre disc (fT3|) is negative for x < 0. Figs. [5]^a) and (b) display level 
curves and the surface plot of the imaginary part of the dimensionless surface 
density with shift parameter a/6 = 0.5. 

Now the combination of the imaginary part of a complex-shifted system 
with another potential-density pair may generate well-behaved and interest- 
ing new configurations. To illustrate our point, we take the superposition of 
the real and the imaginary part of the complex-shifted Kuzmin disc ([6]) 

=K(c7e)+/39(cTe), (43) 

where /3 is a dimensionless parameter which we assume positive, and an 
analogous expression holds for the potential. We have chosen this example 
for its simplicity, but in principle the superposition does not need to involve 



15 




(b) 



lm(a ) 
<= 




Figure 5: (a) Isodensity curves and (b) surface plot of the surface density 
Q{ac)/{M/b'^) equation ([T3|) with shift parameter a/b = 0.5. 



the real and imaginary parts of the same shifted system. Using equations 
(HOll-dlSl), we obtain 



GM (^x + Vx^ + 4a2x2 + 2(3ax 



sup. 



\/2y/x^~T^a?x^y X + a/x^ + 4a2x2 
Mb 



(44) 



a 



sup. 



'2^ - \/^2 + 4a2^2^ + 2/3ax (^2e + Te^Ti^ 



^ + a/^2 _^ 4q_23,2 



(45) 



Isodensity curves of (j45p are shown in Figs. E^a) and (b) for shift pa- 
rameter a/b = 0.5, /3 = 0.5 in Fig. M.a) and /? = 0.8 in Fig. El^b). The 
mass distribution is no more symmetric with respect to the x = axis. For 
X > 0, the imaginary part of the shifted density adds mass to the real part, 
and when x < there is a subtraction of mass. The maximum of density 
is displaced from the centre along the x-axis; in the examples shown in Fig. 
[6)^a) the maximum is located at x/b ^ 0.09 and in Fig. [Q^b) at x/b ^ 0.13. 
For the values of /3 shown in Figs. El^a) and (b), the surface density is non- 
negative everywhere, but one can expect that if /? becomes too large regions 
with negative densities should appear. In fact, we find the following inter- 
vals of P, depending on the shift parameter, where the superposed density 
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Figure 6: Isodensity curves of the surface density OguyJiMl^) equation 
(p5]) for shift parameter ajh = 0.5 and (a) P = 0.5, (b) P = 0.8. 



is non-negative: for a/b = 0.25, < /? < 2.5; for a/b = 0.5, < /? < 1.0; 
and for a/b = 0.75, < /3 < 0.3. 

It is also possible to "inflate" the flat potential-density pair (j44p and 
([15]) by using the same procedure described in Section [5l One obtains the 
following result: 



2Pax 



sup. 



Psup. 



(46) 



4V27r (C^ + 4o2x2)3/2 + TC^TI^ 

^2C - + 4a2x2) + 2pax (2C + ^/C + 4a2x2 
M (1 - hi) 



{-(CV4aV) 



C + + 4a2x2 ) ( 2C - VC^ + 4a2x2 ) + 2(3ax ( 2C + ^/C^^Icfl^ 



Ay/2TT {C + 4a2x2)S/2 ^ y/(2 + 4a2^2 



+^{b + hf 3C^-2Cx/C2 + 4^-4a' 



,2 2 



+2/3aa;(^3C^ + 2C\/C^T4^2^- 40^x2) | 



(47) 



where, as before, ^ = 3:^^ + 2/^ ~ + + h{zy\^ . As an example we take 
h{z) = \J + c^. Figs. [Ha)-(c) show isodensity curves of the mass density 
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-4 -2 2 4 
x/b 



Figure 7: Isodensity curves in the three orthogonal coordinate planes of 
the surface density PsupJ {M /h^) equation ([^7]) for parameters a/h = 0.5, 
(3 = 0.8 and a Miyamoto-Nagai transformation with parameter c/b = 0.3. 



Psup./{M/b^) on the three orthogonal coordinate planes for parameters a/b = 
0.5, /3 = 0.8 and c/b = 0.3. In Figs. El^b) and (c), the values of the level 
curves are the same as indicated in Fig. [7l[|a). Again the lack of symmetry 
with respect to the x = axis is noted. 

7 Discussion 

The method of complexification, also known as complex-shift method, was 
used to generate a family of non-axisymmetric flat discs using the Kuzmin- 
Toomre family of discs as the parent system. The discs present non-negative 
surface densities for restricted ranges of the shift parameter. Further we 
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showed how these discs can be superposed to yield non-axisymmetric struc- 
tures with vanishing surface density at the origin that can be considered 
as complex-shifted flat rings. We also "inflated" the aforementioned flat 
potential-density pairs by a general transformation and presented some ex- 
amples by using a Miyamoto-Nagai transformation and a class of polynomi- 
als. These resulted in triaxial potential-density pairs. We also showed that 
the imaginary part of a complex-shifted system, which has non-physical 
characteristics, can be combined with other pairs to generate well-behaved 
systems that are not symmetrical with respect to a coordinate axis. 

The analytical potential-density pairs presented in this work can all be 
expressed in terms of elementary functions, and yet show non-trivial fea- 
tures. We hope they may be useful as more realistic models in galactic 
dynamics. 
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